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Abstract. Heat-induced mobility of nucleosomes along DNA is an experimentally well-studied phe- 
nomenon. A recent experiment shows that the repositioning is modified in the presence of minor-groove 
binding DNA ligands. We present here a stochastic three-state model for the diffusion of a nucleosome 
along DNA in the presence of such ligands. It allows us to describe the dynamics and the steady state of 
such a motion analytically. The analytical results are in excellent agreement with numerical simulations of 
this stochastic process. With this model, we study the response of a nucleosome to an external force and 
how it is affected by the presence of ligands. 
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Eucaryotic DNA is packaged inside the nucleus by be- 
ing wrapped onto millions of protein cylinders. Each cylin- 
der is an octamer of eight histone proteins and is associ- 
ated with a 147 basepairs (bp) long stretch of DNA [T]. 
The resulting complexes, the nucleosomes, are connected 
via stretches of linker DNA; typical linker lengths range 
from 12 to 70 bp [2\. From the crystal structure [1^ one 
knows that the DNA is bound to the octamer at 14 bind- 
ing sites that correspond to the places where the minor 
groove of the DNA faces the octamers defining the bind- 
ing path to be a left-handed superhelix of one and three 
quarter turns. 

With around three quarters of its DNA being tightly 
bound to the octamers the DNA binding proteins face the 
challenge that their target sites might be masked if they 
happen to be occupied by a nucleosome. One knows of 
two possible pathways of how such wrapped DNA por- 
tions become accessible - at least temporarily. One possi- 
bility is the spontaneous unwrapping of nucleosomal DNA 
from the octamer [31IH[5] that provides DNA binding pro- 
teins a window of opportunity to bind to its target. An- 
other possibility is that the octamer moves as a whole 
along the DNA [g l Tf l lS l ig i fTU l ITT] . thereby releasing previ- 
ously wrapped portions. This so-called nucleosome sliding 
is what we study in the current paper. 

A typical experimental system to investigate nucleo- 
some repositioning consists of a DNA template slightly 
longer than the wrapped portion (e.g. 207 bp in Ref. [7]) 
that is complexed with one octamer. The nucleosome po- 



sition can be inferred from the electrophoretic mobility of 
the complex in a gel. It is found that nucleosome sliding 
is a slow process and that it takes a nucleosome around 
an hour to reposition completely on such a short DNA 
fragment. Another important observation is that the new 
positions are all multiple of 10 bp (the DNA helical repeat) 
apart from the starting position. 

Concerning nucleosome sliding there are currently two 
possible mechanisms discussed that could explain those 
experiments [9lll0j. Both mechanisms have in common 
that they rely on defects that are thermally injected into 
the wrapped DNA and that traverse the nucleosome thereby 
causing its displacement. The reason for assuming defects 
as the cause of repositioning rather than the sliding of the 
octamer as a whole is that the latter mechanism would re- 
quire the simultaneous detachment of all 14 binding sites 
which is too costly (around 85 ksT). The two kind of de- 
fects are 10 bp loop defects [l2lfT3] and 1 bp twist defects 
[I4|[T5]. A 10 bp loop defect is a bulge that carries an ex- 
tra length of 10 bp causing redistribution events of that 
step length. These preserve the rotational orientation of 
the nucleosome. This fact as well as the predicted value 
for the mobility seems to agree with experiments [13] . 

The second class of defects, the 1 bp twist defects, 
carry a missing or an extra bp. To accommodate such a 
defect between two nucleosomal binding sites the DNA 
needs to be stretched (or compressed) and twisted (hence 
the name) . A nucleosome mobilized by twist defects moves 
via 1 bp jumps. Since the octamer is always bound to 
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the minor groove of the DNA, the nucleosome performs 
a corkscrew motion around the DNA. Alternatively one 
can say that the DNA acts as a molecular corkscrew. 
Since twist defects are much cheaper than loop defects 
(« gfcsT [H] vs. « 23kBT [13]), twist defects are expected 
to make nucleosomes much more mobile than observed 
in experiments. However, repositioning experiments are 
always done in the presence of so-called nucleosome po- 
sitioning experiments that are now known to be wide- 
spread in eukaryotic genomes \13^. Such positioning se- 
quences (like the sea urchin's 5S rDNA sequence of Ref. 
[7j) make use of the fact that certain bp sequences induce 
an anisotropic bendability of the DNA. A nucleosome dif- 
fusing via twist defects feels the resulting energy landscape 
since in order to move by 10 bp along the DNA it needs 
to perform a full 360° rotation around the DNA (unlike 
a nucleosome that moves via 10 bp bulges). In the end, 
both mechanisms are predicted to appear very similar in 
the experiments mentioned above |10| . 

There is an experiment [T7j that hint at twist defects 
being responsible for nucleosome sliding. This experiment 
is performed on a 216 bp DNA that contains the sea urchin 
5S rDNA sequence in the presence of minor groove bind- 
ing pyrrole imidazole polyamides, synthetic ligands that 
can be designed to bind to short specific DNA sequences. 
It was found that the nucleosome mobility is dramatically 
reduced when such ligands are added. This might reflect 
the fact that a DNA corkscrew motion is sterically forbid- 
den once a ligand is bound |15| . 

Finally, the nucleosomal mobility might also be im- 
portant when a transcribing RNA polymerase encounters 
a nucleosome. Experiments with short DNA fragments 
carrying single nucleosomes show that in such a setup a 
transcribing RNA polymerase (e.g. bacteriophage poly- 
merase from T7 [T7] and SP6 [H]) can transcribe the 
whole fragment, even though it is partially occupied by 
a nucleosome. An interpretation of how the polymerase 
negotiates with the nucleosome is tricky since there are 
at least two possible explanations. The simpler explana- 
tion is that the polymerase crosses the nucleosome in a 
loop [TSlfTO] . However, the alternative explanation [15^^ is 
that the polymerase pushes the nucleosome in front of it, 
pushing it off the template; before the octamer falls off it 
rebinds at the other DNA terminus. Interestingly, in the 
presence of ligands the polymerase stalls [17] , pointing to- 
wards the second mechanism. Note that the former mech- 
anism would also allow a polymerase to negotiate with an 
array of nucleosomes, the second does not ("traffic jam"). 
In fact, experiments [20] show that RNA polymerase can 
only transcribe through an array and leave it intact if a 
nuclear cell extract is present, otherwise the nucleosomes 
are stripped off the DNA. 

In this paper, we study the influence of ligands on the 
mobility of a nucleosome along a long sequence of DNA 
that is either assumed to be uniform ("random DNA") 
or periodic ( "nucleosome positioning sequence" ) . First we 
solve a three-state model inspired from Ref. [15] using 
a method described in [^[^ . There are some data ex- 
tracted from experiments in which the mobility of the nu- 
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Binding site occluded 
Binding site open 
Ligand bound 




AG = e/fcfiT 



Fig. 1. Schematic view of the three-state model for sliding of a 
nucleosome on a periodic DNA sequence of height AG = eksT 
and period 21. States 0, 1, and 2 denote states where the ligand 
is bound, unbound with its binding site open and unbound with 
its binding site closed, respectively, a;^^ ' denote the rates of 
going from state i to state j to the right (left). 



cleosome was probed using electrophoresis methods f l7) . 
The results of the model are then compared with these ex- 
isting data for short DNA. Finally, we study the effect of 
an external force on the sliding of the nucleosome in this 
model. In a very simplifled way this force can represent 
the action of an RNA polymerase or other DNA binding 
proteins on a nucleosome [15] . 



1 A stochastic model for nucleosome sliding 

In the absence of ligands, the sliding of the nucleosome 
along its wrapped DNA can be described by a hopping 
model with a single state or with two-states depending 
on the DNA sequence. The effect of the sequence can be 
modelled as shown in Fig. [T] (without considering state 

in this figure). State 1 represents the preferred binding 
sites for the DNA-histone complex (the minimum of the 
sequence dependent potential), while state 2 may be the 
high energy state of this potential. For uniform sequences 
the height of the potential barrier is zero (the energies 
of the two states are equal) that corresponds to a single- 
state model. For a positioning sequence of DNA, states 

1 and 2 are separated by 5 bp, i.e., half of the period of 
the sequence denoted as I in the figure. Similar two-state 
models have been used to describe the motion of a linear 
molecular motor on its track (such as a kinesin walking 
on a microtubule) [^T[|22j . 

The presence of ligands affects the interaction of the 
nucleosome with the DNA, and we describe this as an ad- 
ditional state that branches off the state 1 or 2. When 
the ligand binds to the DNA, we assume that no sliding 
can happen, as confirmed by experiments [17]. So the nu- 
cleosome waits until the ligand detaches from the DNA, 
and this will reduce the mobility of the nucleosome. The 
ligand may bind on the DNA either in state 1 and 2 de- 
pending on the location of its substrate that is a short 
DNA sequence. In this model we only consider the case 
where the ligand can bind to the DNA when the nucleo- 
some is in state 1. This is equivalent to AG > in Fig.[T] 
In state 2 (5 bp, half the DNA pitch, away from state 1) 
the binding site is then inaccessible since it faces the oc- 
tamer surface. Similar periodic sequential kinetic models 
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with branching, jumping or deaths have been studied in 
[23] , Especially in Ref. [15] the three state model in Fig. 
[1] has been examined in the same context as the current 
paper for special cases of the transition rates. 

To model the dynamics of the histone, consider the 
one-dimensional lattice of Fig. [TJ The nucleosome can hop 
to neighboring sites on this lattice with some specific rates. 
Assume that I is the distance between sites 1 and 2 (which 
is half the period), and let pi{n,t) be the probability for 
the nucleosome to be in the states i = 0,1,2, at posi- 
tion X = nl and at time t. The pi{n,t) satisfy the master 
equation 

dtPoin,t) = ujioPi{n,t) - ujQipa{n,t) (1) 

dtPi{n,t) = p2{n - l,t) +uj2iP2{n + l,t) 

+ LL>oipo{n,t) - (w+ +UJ^2 +(^io)pi{n,t) 

(2) 

dtP2{n,t) = cjj*2Pi(" - 1,^) +^i2Pi{n' + 1)^) 

- {uj+^+uj2i)P2{n,t), (3) 

where cj^j i^Tj) represents the rate of transition from state 
i to neighboring state j on the right (left); and wio (^oi) 
represents the binding (unbinding) rate of a ligand. 

Let us introduce a vector F{X,t), the components of 
which are generating functions of the position for each 
state i. These components are Fi{X, t) — ^ e~^"'''pi{n, t). 

n 

The use of generating function reduces the space of con- 
figurations that is infinite, to only 3 states due to the 
periodicity of the system. It makes the calculations much 
more tractable. The master equation now becomes: 



and 



dtF{X,t)^M{X)F{X,t), 



(4) 



^ _ 1 rf^cr^ 



2 dX^ 



(8) 



A=0 



One can expand o'm(A) near A = as crm(A) = —vX + 
DX"^ + 0{X^). Using this expansion in the eigenvalue equa- 
tion, det [AI — (Tm(A)/3] = 0, the velocity and the diffusion 
constant are derived as 



_ ^ 2 (^]^2^21 ~ ^12^21) ^ 

S+^iu;ti+u;2i) ' 



(9) 



D = 2V 



42^tl + ^^12^^21) K + 8tJ]*2^12^^21^21 + 



with 



S = ^++^^12+^^1 +^21. 



-r, +2 _2 i2 _2 



(10) 



(11) 



+ 2 (^^1^2^12 + ^12^21 + ^21^21 + ^21^12) . (12) 



(w+ -I- W21) 



— 2 (lJJ^2^21 ^ ^12^21) 



21 



1 - 



^21) 



(13) 



with 



It is worth to mention that in the ligand-free case the rate 
of going from site 1 to is zero, loiq = 0, and the model 
becomes equivalent to the two-state model that has been 
discussed in j21,22j . 



M(A) 



-^01 



-'12 



'12 



- uJiQ UO21 e ^' + LU21 e^' 



Wj*2 e~-^' -I- uj 



12 ^ 



-OJ. 



21 



'21 



The solution of Eq. ([4]) is 

F(A,i) = e*^(^)*F(A,0). 



(5) 



(6) 



After calculating the eigenvalues of the matrix A/, one can 
see that at long time t +00 only the largest eigenvalue 
of M that is denoted by (7™ (A), contributes to -Fi(A, t) = 

i 

^g-AnZ\^ ~ g'^™(A)t_ jv^g^g i^T^^^ ^j^g normalization condition 
for the probability implies that (7^(0) = 0. The eigenvalue 
Cm (A) contains all the long time dynamical properties of 
the system, such as the velocity and the diffusion constant 
[2T1[^. since 

v = ^ , (7) 



dX 



2 Kinetic rates in the absence of force 

We assume that the nucleosome sliding in the absence of 
any external force is a passive process, which implies that 
there is no preference between left and right direction. 



Thus uj^2 — ^12 ~ "^12 UJ21 
introduce k as the ratio of rates 

(^12 
K = — e 

^21 



UJ- 



21 



0^2 1. Let us 



(14) 



A=0 



where e = AG/ksT is the energy difference between the 
two states 1 and 2, and the above equality expresses the 
detailed balance condition. 

We consider the binding chemical reaction of the lig- 
and: L + S ^ LS, where S is the substrate, i.e. the nucleo- 
somal DNA, L is the ligand and LS is the ligand-substrate 
complex. From kinetics theory of first order chemical re- 
actions, one can write luiq = k+[L][S] and wqi = k-[LS], 
such that the equilibrium constant of the chemical reac- 
tion is Keq = k+/k- = [LS]eq/ {[L]eq[S]eq), in tcrms of 
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equilibrium concentration of ligand [L\eq, substrate [S]eq 
and complex [LS'Jgg. Therefore the ratio 

^c^^ [L][S] _ [LSU [L][S] 
^~ iooi [LS] [LU[SU [LS] ' ^ ^ 

quantifies the deviation away from equilibrium. In general, 
the substrate is in excess so that [S] and [S]eq are both 
large, and also [S] ~ [S]eq- Furthermore, if one assumes 
that [LS] ~ [LS]eq, then 77 ~ [L]/[L]eq. 

From Eq. one finds that w = as expected and 
the diffusion constant is: 



D 



2 UJ12 f 



1 



(16) 



where 0^12 can be calculated by Kramers rate theory in 
the limit e ^ 1. 

The sequence dependent potential can be modelled as 
a periodic function j3U{x) = | cos (^) with /3 = l/ksT 
[14]. The time needed for the nucleosome to go from one 
minimum of this potential to the neighboring maximum, 
T, can be determined using the Kramers rate theory : 



1 



Do 



2TTkBT 



^\U"il)U"{2)\ = -\^Do 



4P 



(17) 



where Do is the diffusion constant of the nucleosome with- 
out any potential barrier or ligand, 1 and 2 refer to the 
minimum and the maximum of the potential U{x). From 
this we find uji2 for such a strong positioning sequence: 



1 TTe _ 
Yr = 8l^^° ' 



for e > 0, 



(18) 



i being the probability to go through the 



with the factor 

barrier 2 from either direction [25 • It is worth to mention 
that in the limit e ^ 1 with (e > 0, 77 > 0) when there is a 
sequence dependent potential for the nucleosome, one can 
find that: 

In the case of random sequence of DNA, i.e in the absence 
of a sequence dependent potential, e = and then uji2 is 
simply equal to 

Do 



W12 



and the diffusion constant can be found as 



D 



2Do 
2 + 77' 



(20) 



(21) 



Putting in numbers, using realistic parameter values [14j 
Do — 600 bp^/s, e — 9,1 = 5bp,ri = in the absence of lig- 
and, and rj ~ 100 in the presence of ligands our model pre- 
dicts the time needed for a nucleosome to diffuse on a 70 
bp DNA, to be 78 minutes without ligands and 131 hours 
with the ligands which is consistent with the experimental 
observations [T7]. For a random sequence (e = 0), in the 
presence (77 ~ 100) and the absence of ligands (77 ~ 0), the 



characteristic time for 70 bp diffusion are 3.5 min and 4 
sec, respectively. 

In the absence of ligands and for arbitrary sign of e, 
Eq. becomes 



D^^. Do 



(22) 



The above expression that has been derived with a discrete 
stochastic model, can be also derived using a continuous 
description in the limit e ^ 1. This can be done by con- 
sidering a particle diffusing in the periodic potential U (x) , 
for which the diffusion coefficient is |!24 



D = Do{e^"'^''Y\e-P''^^Y\ 



(23) 



Indeed (e'^'^^^)) is a Bessel function, the asymptotic form 
of which is 2e''/'^ / for e > 1. From this Eq. ^ is 
recovered using Eq. ([23]) . 



2.1 Transcription-induced sliding 

Up to now we considered thermally induced, undirected 
nucleosome sliding. Here we discuss the case when a force 
is applied to the nucleosome. As discussed in the introduc- 
tion this situation might occur when an RNA polymerase 
encounters a nucleosome during transcription. A more mi- 
croscopic model of such an encounter is presented in Ref. 
[25]. Note that we will assume that the polymerase does 
not unpeel the nucleosome, a case considered recently by 
T. Chou |26j . Rapid progress in the field of micromanip- 
ulation experiments let us expect that there will be also 
soon data available where forces are directly applied to 
nucleosomes. 

A force F exerted on the nucleosome introduces a bias 
in the transition rates: 



9tf 



9r/ 



(24) 
(25) 
(26) 

.2l-we-^ (27) 

where Of are the load distribution factors [23], and / = 
Fl/kBT. Using detailed balance condition, one has df + 
^1 +0t +0^ ^2. 

Putting these rates in Eq. ([9]), the velocity of the nu- 
cleosome can be found as 

V — 2ujI e 7 r 7 r- 

e-' (e»^/ -I- e-»r/ j -f (1 -I- 77) (e<^tf + e~«2"/j 

(28) 

The mobility of nucleosome is defined as 



1 dv 



keT df 



which gives 



1 



2ajK 



-I' 



ksT 1 + 77 + K 



(29) 
(30) 
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F(pN) 

Fig. 2. The velocity of the nucleosome versus the external 
force exerted on the DNA for (a) a random sequence of DNA 
(e = 0) and (b) a positioning sequence with e = 9 in two cases: 
the diamonds are simulation data for -q = while the black 
solid line is plotted using theory, Eq. (|28|l . the triangles are 
simulation data for rj = 100 while the brown dashed line is 
plotted using theory, again Eq. (f28|l . 




100 



Fig. 3. Velocity v of the nucleosome versus rj in the presence 
of an external force F = 10 pN on (a) a random basepair se- 
quence (e — 0) and (b) a positioning element (e = 9). The red 
diamonds are the simulation results while the line is plotted 
using the analytical approach, Eq. 1281 



Comparing this expression with Eq. ()16p one finds the 
Einstein relation = D/^ksT) verified. 



3 Results 

There are three physical quantities that affect the behavior 
of the system: the external force F, the sequence depen- 
dent part of the potential measured by e and the ligand 
concentration [L] that enters into rj through rj ~ [L]/[L]eq- 
The physical behavior of the system is characterized by the 
velocity of the nucleosome repositioning along the DNA 
and its diffusive behavior. In this section, the results of the 
analytical approach described in the previous sections and 
of a computer simulation that is discussed in Appendix A, 
are presented. 



We first consider the effect of an external force on the 
velocity of the nucleosome repositioning along the DNA. 
In Fig. [2] we plot the nucleosome velocity v versus the 
applied force F for the two limiting cases 77 = and 77 = 
100 of the experiments [T7], on both random DNA and on 
a positioning sequence. As expected, the velocity of the 
nucleosome increases with the external force and there is 
a good agreement between the simulation results and the 
analytical approach. At zero force, the nucleosome shows 
purely diffusive behavior and there is no net velocity, cf. 
Fig.[2]at F = 0. As soon as a force is applied, there is a bias 
in the transition rates and the nucleosome attains a drift 
velocity in the direction of the applied force. A positioning 
sequence of DNA leads to an effective potential barrier on 
the corkscrew path of the nucleosome [M] leading to a 
drift that is significantly smaller than on random DNA. 
The aforementioned behaviors are seen in Fig. [2l 
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Fig. 4. Diffusion constant D of the nucleosome versus external 
force _F on a positioning sequence (e = 9). The solid black 
line (red diamonds) and the brown dashed line (light brown 
triangles) are the theoretical (simulation) results for = and 
t] = 100. 



Next we study how the hgands influence the shding 
velocity. In Fig. [2] we see that for a given finite force the 
nucleosome velocity for the case r\ = 100 is much smaller 
than in the absence of ligands, 77 = 0. The effect of the lig- 
and concentration on the drift velocity of the nucleosome 
for a typical external force, F = 10 pN, is shown in Fig.[3l 
As expected ligands block corkscrew sliding, lowering the 
overall drift velocity. For typical experimental numbers 
|17j already a small concentration of ligands significantly 
lowers the drift, cf. Fig. [31 

Another parameter that provides information about 
the system is the diffusion constant of the nucleosome, D. 
The behavior of D versus both the external force and the 
ligand concentration has been checked on for random DNA 
and a positioning sequence. In Fig. [4] it can be seen that D 
increases with F. Using a simple two state model can help 
us to understand the behavior of the diffusion constant in 
terms of force F. This simple model has been explained 
in Appendix B. We see that in the presence of external 
force, the diffusion constant of the particle becomes larger 
when the force is increased. 

We present the behavior of the diffusion constant ver- 
sus 7], for two different forces, F — OpN and F = 10 pN, 
in Figs. [5] and [HI Naively one would expect that at a fixed 
external force the diffusion constant decreases with 77 since 
an increase of the ligand concentration leads to a higher 
probability to have a ligand bound that then suppresses 
diffusion. For zero force the diffusion constant does in- 
deed follow this expectation, cf. Figs. [5]Ja) and[SJa). In- 
terestingly, in the presence of a nonzero external force the 
behavior of the diffusion constant versus rj differs dramat- 
ically from this expectation. For rj < 1, D increases with 
ry, and then decreases as rj goes to infinity (Figs. [5jb) and 
inijb)). For random (positioning) DNA the maximal value 
of the diffusion constant is five (two) orders of magnitude 
larger than the value in the absence of ligands. 




Fig. 5. Nucleosome diffusion constant D versus rj for random 
DNA (e = 0) for two different forces: (a) F = OpN and (b) 
F — 10 pN (the inset gives a zoomed view showing the non- 
monotonous behavior for small 77- values). The red diamonds 
are the simulation results while the lines use the analytical 
approach, Eq. 1101 



For 77 ^ 1 the diffusion constant in 3-state model be- 
haves as the 2-state model with a correction like D3_statc = 
-D2-statc + 0(77). For large values of 77, as can be seen in 
Eq. (fTOl) . the diffusion constant changes as D oc I/77 and 
the diffusion constant decreases as 77 increases. The 77 at 
which D attains its maximum 77* can be calculated from 
Eq. ([To]) and is shown in Fig. [7] as a function of F for two 
cases of random and positioning DNA. As can be seen in 
this figure, at not very small forces, 77* is equal to 1 for ran- 
dom DNA and 0.5 for the positioning sequence. Since at 
large forces the positive rates dominate, the diffusion 
constant, Eq. p^ . simplifies to 



D = 2f- 



77(w+ 



12'^2i; 



^01 [^^2liV + 1) + 



4]^ 



(31) 
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Fig. 6. The diffusion constant of the nucleosome versus rj in 
the case of a positioning sequence with e = 9 for two differ- 
ent forces: (a) F ^OpN and (b) F = lOpiV. The inset is the 
zoomed-in plot showing the behavior of the diffusion constant 
for small values of rj. The red diamonds are the simulation re- 
sults while the lines are plotted using the analytical approach. 



Setting the derivative of D with respect to rj equal to zero, 
one obtains 



2lu+ 

1 + e-^+^'^i~'^2)f 



(32) 



where we have used Eqs. (IM])-([27|). For the case df" = 1/2 
the force dependent term drops out and we find rj* = 
(1 -|- e~'^)/2 and especially 77* = 1 for random DNA. We 
shall come back to this point in the discussion section. In 
Fig. [71 the behavior of rj* is plotted versus the external 
force in the case of 9f = 1/2. 

Let us now discuss the behavior of the diffusion con- 
stant versus rj for small values of 77, i.e. for 77 ^ 1. Through 
an expansion of the exact expression we obtain 
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Fig. 7. rj* versus F for the case Of" = 1/2 for random DNA 
(dashed line) and a nucleosome positioning element (solid line). 
The lines are drawn using Eq. 1321 



where Aq and Ai are functions of / that are given in Ap- 
pendix C. It is convenient to expand Ai{f) as a function 
of/: 



dD 
drj 



aiPo+Pif) + 0{f), (34) 



77=0 



where a, (3o, and Pi are functions defined in Appendix C. 
Since a > and /3o < 0, we have dD/dr] < for both ran- 
dom and positioning sequence of DNA at / = that can 
indeed be seen in the plots of diffusion constant D versus 
?] shown in Figs. (Ha) and [51^ a). When /?i becomes posi- 
tive, there is a threshold in force that is fr — \/—fh/Pi 
such that for f > fr the derivative of D with respect to 77 
is positive for small 77-values. From the expression of the 
coefficients derived in Appendix C, the value of /t is the 
following function of e, uj and ojqi: 



V2{1 



-3 + e-^ (5 - 6e-^) -I- 4e-« (1 + e~^) ■ 



1/2 



(35) 



D^Ao{f)+Aiif)r, + 0{v^), 



(33) 



Putting in numbers, ujqi — 0.001 s^^ and uj — Dq/P ~ 
24 s""'^, we find /t — 0.48 that is equivalent to F ~ 
0.016 pN for a random sequence of DNA. In the case 
of positioning sequence of DNA the rate of uj changes to 
w = TreDo/{8P) ~ 85 s'^ and we find fr 0.23 that is 
equivalent to F ~ 0.55 pN. In Fig.[7]the behavior of 77* as 
a function of F is shown for both random and positioning 
sequence of DNA. We note that the force value for which 
77* = corresponds to the threshold force given by Eq. 

How can the surprising non-monotonous behavior of 
D as a function of 77 of a nucleosome driven by the ap- 
plication of a force be explained, especially the strongly 
enhanced fluctuations around 77* with D{rf) ^ _D(77 = 0)? 
Obviously the fluctuations of position of this driven nu- 
cleosome in the presence of ligands are of different origin 
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than the ones in the absence of ligands. For sufficiently 
large forces the nucleosome mostly steps in the direction 
of the force or - if the nucleosome is in state 1 - a ligand 
might bind. The latter event stops the drifting nucleo- 
some for a while and is thus a source of fluctuations of 
completely different origin, because this introduces some 
waiting time before hopping from states 1. The higher the 
concentration of ligands, i.e. the higher the value of 77, the 
more often these events occurs, increasing their contribu- 
tion to the overall fluctuation of the nucleosome, which are 
measured by the diffusion coefficient. This is the case up to 
a critical value of ry, named 77*, which is force dependent. 
Further addition of ligands populates the state so much 
that the nucleosomes gets frequently stuck, decreasing its 
diffusion constant. 



4 Discussion 



In our model the external force changes the local rates 
from sites 1 to 2 and vise versa. This is because the force 
produces internal stress on the nucleosomal DNA that in- 
troduces a bias in the dynamics. This assumption can be 
verified by considering the microscopic details of the in- 
teraction of the nucleosome with the DNA that will be 
presented in a forthcoming publication [25j . 

In the above treatment we assumed that both binding 
and unbinding rates of ligands to the DNA do not depend 
on the external force. The typical length of a ligand site is 
about 6-7 bps and the DNA length between two adjacent 
binding sites is 10 bps. Even at the highest forces consid- 
ered here (30 pN) the probability of having a defect in a 
ligand binding region is close to zero [25] so that we expect 
that its influence on the ligand rates can be neglected. 

It is also important to point out that the load distribu- 
tion factors of have an effect on the diffusion constant. We 
assumed that all Of are equal to i. By changing the val- 
ues of of the overall behavior of all plots does not change, 
although the precise values of the diffusion coefficient and 
even the curvature of the plots can be affected. For in- 
stance, the value of ry* at which D is maximized depends 
on the values of Of as can be seen from Eq. ([32]) . The mi- 
croscopic details of the interaction between the DNA and 
the nucleosome determine the values of the distribution 
factors. For a random sequence there is no reason to have 
different rates for backward and forward steps of the nucle- 
osome along the DNA. Therefore, one finds that for large 
values of the external force, 77* converges to 1. For a DNA 
positioning sequence, the values of 0^ could in principle 
depend on the sequence. The exact value of these coeffi- 
cients can only be determined from experimental data or 
from a more detailed modelling of the transition state. We 
have arbitrarily chosen them to be 1/2 for the plots. Note 
that if of < of then ri*\p^oo = 0.5, the same value as in 
our case 9f = Of , while for Of > Of , rj* increases as F is 
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6 Appendix 

A The algorithm of the simulation 

The 3-statc model presented in this paper is simulated 
using a "Random Selection Method" [27j. It is defined in 
terms of the transition rates tof^ that give the probabili- 
ties per unit time for going from state i to state j in the 
plus/minus direction. If the system is at time t in the state 
i, a transition to the neighboring state j happens at time 
t + At with the finite probability — At ujfj . 

For each step, a random number < < 1 is drawn. 
Depending on its value and the state of the system a de- 
cision is taken: 

state: 0^1 if < C < Poi, 
state : ^ otherwise 



state :l->0 ifO<C< Pio, 

state : 1 ^ 2+ if Pio < ^ < Pio + P^2: 

state : 1 ^ 2" if Pm + P]*2 < ? < Pio + P^- + Pr2 

state : 1 — > 1 otherwise 



state : 2 ^ 1+ if < ^ < PJl, 

state : 2 ^ 1" if P^i < ^ < PJl + P21, 

state : 2 ^ 2 otherwise. 

For the next time step, from t + At to t + 2At, the pro- 
cedure is repeated again. The time step At is chosen small 
enough, such that for each step the condition J2 ^fj^'^ < 1 
is satisfied, where the sum is taken over the probabilities 
of all possible transitions from state i. 

This algorithm is similar to the one of Gillespie [28], 
except for the fact that we use here a constant time step 
whereas for the Gillespie algorithm the time step is a 
random variable. Both algorithms converge to the same 
steady state albeit after different times as we also checked 
for our model. The steady state probabilities for the three 
states are obtained by setting the time derivatives of the 
probabilities in Eqs. ([J)-© to zero: 



increased and one has r]*\p- 



00. 



Pq 
Pi 

P2 



1 - 


h 7/ + e ^ 




1 


1 - 


1- 77 -f e^' 







1 -I- 77 -I- e 



(36) 
(37) 
(38) 
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Fig. 8. The two state problem used to explain the behavior 
of the diffusion constant versus the applied force. The external 
force, F, changes the jumping rates of the particle, which are 
shown by g{F) and r{F). 



a result that has been previously obtained in Ref. [15'. 
We let the simulation run for a long time (from to to 
tN with ti — iAt, ^ 1) to be sure that the system 
has reached equilibrium. Then averaged over M ensem- 
bles with M ^ 1, the mean velocity and the diffusion 
constant is determined by 



Xrih)) 



NAt X M 



(39) 



and D is determined as the slope of the plot of X^rli X'^ /M- 

(12^1 Xr/Mf versus 2 Y.r=l ^r/M. 

The time steps used for the simulation are less than 
0.001, depending on the simulated case. The time goes to 
10^ s, and the number of ensembles are M = 2000. The 
used parameters for the simulation are 6*^^ = 1/2, i — 1,2 
and Lo is determined from to Eq. (fT8|) and Eq. (pO| for the 
positioning and random DNA sequence respectively. Also 
using the experimental data, the typical time needed for 
a ligand to unbound from the DNA is some minutes and 
woi = 0.001 s-i. 



B The behavior of the difFusion constant 
versus the external force in a simple two 
state model 



increases and one decreases. Here, we have assumed that 
the force pushes the system to the right, so g{F) increases 
with F while r{F) decreases with F. Then a simple cal- 
culation leads to: 



d_ 

dt 
d_ 

dt 



(n^) - 2 [g{F) - r{F)] (n) + giF) + r{F), (42) 



{n)^^2[g{F)-r{F)] (n). 



where we have used: 



dt 
d 



d 



^-(A)^ = 2(A)-(A„), 



and 



(43) 

(44) 
(45) 

(46) 



Using Eqs. (|42| and (|43|) . one finds the diffusion con- 
stant as following: 



D 



g{F) + r{F) 



(47) 



If the rates in the absence of the external force are de- 
noted by w, the external force, F , changes the jumping 
rates of the particle to g{F) = oj cxp[-\-Fa/ {iksT)] and 
r{F) = LO exp [-Fa/ {2kBT)\. Using Eq. ^ and the men- 
tioned rates in the presence of external force, the diffusion 
constant is found as 



^+Fa/2kBT I -Fa/2kBT 



(48) 



This explains the behavior of our results for the diffusion 
constant mentioned in the text. 



To check the effect of force on the diffusion constant, let 
us assume a simple two state model in which the external 
force changes the rates as shown in Fig. [8] The position 
of the particle in the mentioned lattice model is denoted 
by n. By definition, the diffusion constant is written as 



(40) 



where (A) denotes the average of quantity A that is given 
by (A) = J2n ^nPn- Pn IS the probabihty for the particle 
to be in the position n. 

The master equation governing this system can be 
written as 



dPn 

dt 



r{F)F,+, + g(i^)P„-i - [riF) + g{F)] P„, (41) 



where the force is denoted by F, a is the jump length, and 
g{F) and r{F) are the force-dependent rates for going to 
the right and left, respectively. Depending on the direc- 
tion the force is exerted on the system, one of these rates 



C The explicit forms of the auxiliary functions 

In this appendix we give the explicit form of the constants 
used in the Eqs. ((33)) and (p4|) . First we expand D for 



where 



D = 



Bp + B,r^ + Bsf]^ 

[5 + ^27;] 3 

Bq BiS — 3BqB2 



?7 + 0(7?2), (49) 



Bo = [K{uj+2i^+^ + ^^12^21) + 8t^]*2^2l'^12^2l] > 

Bi 2f[2{uj'l2Ujtl +'^r2^2l)('^21 +'^21)5' 

-2(C.+ C.+ - C.r2-2-l)^(l - i^^^l±^)], 

B2 = C^21 + ^21' 

B3 = 2P [(^21 + UJ2lf{ujt2^tl + ^r2^2l)] • 
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From this follow the expansion coefficients in Eq. ([33| : 

(50) 



Mf) = f , 



Mf) = (51) 

Finally we provide here the behavior of Ai{f) for small 
forces. Using Eqs. with 0f = 0.5, the expansion 

of Bi{f) for small forces can be written as 

5 = 2c.(e-^ + l)(l + ^), 

Bo = 8uj^e-'{e-'' + 1)^ 

+ Auj^e-' [2 + 2e-2^ + [e" + if] /^ 

|i = 16c.4e-(e-^ + l) 



16iJ^e-^' ( 1 



24 



21 

Consequently, 



^-2.(1 + ^). 



e-'{-2(l + e-')' 



-3 + 6-'= (5-6e 



4e~Ml + e" 



^01 



(52) 



Using these expansions and Eq. (|34|) . one can write: 

72, ,5 



/3o = -2(1 + e-^)^ 



(53) 
(54) 

/3i = -3 + e~' (5 - 6e-'-) + Aer' (l + e"') — . (55) 
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